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We present a theory of the scaling behavior of the thermodynamic, transport and dynamical 
properties of a three-dimensional metal at an antiferromagnetic (AFM) critical point. We show how 
the critical spin fluctuations at the AFM wavevector q = Q induce energy fluctuations at small q, 
giving rise to a diverging quasiparticle effective mass over the whole Fermi surface. The coupling of 
the fermionic and bosonic degrees of freedom leads to a self-consistent relation for the effective mass, 
which has a strong coupling solution in addition to the well-known weak-coupling spin-density-wave 
solution. We thereby use the recently- introduced concept of critical quasiparticles, employing a 
scale dependent effective mass ratio m*/m and quasiparticle weight factor Z. As a consequence of 
the diverging effective mass the Landau Fermi liquid interaction is found to diverge in all channels 
except the critical one, causing important vertex corrections. The ensuing spin fluctuation spectrum 
obeys c*;/T-scaling. Our results are in good agreement with experimental data on the heavy fermion 
compounds YbRli2Si2 and CeCue-rcAux assuming 3D and 2D spin fluctuations, respectively. 



Introduction. Quantum phase transitions in metallic 
compounds have attracted considerable interest over the 
last two decades. These systems exhibit deviations from 
the standard Fermi liquid model. This "non-Fermi liq- 
uid" behavior is a consequence of the interaction of the 
fermionic (Landau) quasiparticles with bosonic critical 
fluctuations. Early theories of quantum critical behav- 
ior, formulated in the framework of a Ginzburg-Landau- 
Wilson action of the order parameter field <f> [U [5] , found 
that the effective dimension of the corresponding </> 4 -field 
theory is increased to d e f / = d + z where d, z are the 
spatial dimension of the fluctuations and the dynamical 
critical exponent, respectively. In many cases of interest 
d e ff is above the upper critical dimension, so that the 
fluctuations are effectively non-interacting and the the- 
ory is of the Gaussian type (for a review see [3]). While 
this theory is well founded in the case of non-metallic 
systems, for metallic systems the question arises whether 
the fermionic degrees of freedom may be easily integrated 
out. In this paper, we show that it is often not possible to 
reduce the description to a Hertz-Millis model. Rather, 
the interplay of fermionic and bosonic degrees of free- 
dom generates critical behavior of the fermionic quasi- 
particles, acting back on the spectrum of bosonic fluctu- 
ations. Then, a strong coupling regime with respect to 
the fermion-boson coupling is reached, while the interac- 
tion of critical bosonic fluctuations with each other is still 
weak, such that the theory remains in a "renormalized 
Gaussian regime" . 

Experimentally well-studied candidate systems we 
shall focus on are the heavy-fermion compounds 
CeCug-xAuz (CCA) for which we postulate quasi-two- 
dimensional antiferromagnetic (AFM) spin fluctuations 
at the critical point x=0.1 |4j and YbRh 2 Si 2 (YRS), for 
which we assume AFM fluctuations at the critical point 
at magnetic field H = H c of three-dimensional charac- 



ter at temperatures T less than 0.3 K, crossing over into 
three-dimensional ferromagnetic fluctuations at higher T. 
[B] An issue that has hampered progress in developing a 
strong-coupling theory of criticality that involves AFM 
fluctuations has been that they transfer a large momen- 
tum of the order of the ordering wave vector q = Q . 
As a consequence, the self-energy in one-loop order be- 
comes highly anisotropic, being critical only at so-called 
"hot spots" on the Fermi surface that are connected by 
Q. However, the combined exchange of two AFM fluctu- 
ations, which may be viewed as a spin-exchange energy 
fluctuation [7j, may transfer only a small momentum q 
and we argue here that the one-loop order process of such 
energy fluctuations is dominant in providing a renormal- 
ization of the fermionic quasiparticle effective mass that 
is uniform over the Fermi surface. Thus, we consider 
here the simplest case, in which the fermion self energy 
is only weakly momentum dependent. A different way 
by which the effect of the critical AFM fluctuations may 
be distributed all over the Fermi surface is by means of 
impurity scattering [5J [S] ■ 

However, a problem with multiple fluctuation exchange 
is that each additional fluctuation propagator gives rise 
to an additional energy integration and will thus con- 
tribute a small phase space factor, anywhere from lo 2 to 
ui, depending on the critical momentum dependence. As 
we will show below, such factors may be offset by Fermi 
liquid renormalization factors Z~ n >> 1, provided the 
quasi-particle weight factor Z tends to zero, as the exci- 
tation energy uj or the temperature T tends to zero. As 
discussed in our previous work [BJ [9] , even if Z — > in 
the limit u, T — » the quasiparticle picture may still be 
applicable at non-zero u, T, since the quasiparticle width 
r gets renormalized by a factor of Z, which helps to keep 
it smaller than the excitation energy u, a necessary con- 
dition for the existence of a quasiparticle peak in the 
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single-electron spectral function. A careful identification 
of the renormalization of the bosonic (AFM) propagator 
by Z-factors is at the heart of our theory. This includes 
the determination of critical vertex corrections of various 
types. For example, while it is well-known that the spin 
susceptibility at small q receives a correction of the ex- 
ternal vertex of order 1/Z, which follows from a Ward 
identity, we argue that a similar vertex correction ap- 
plies even at the AFM wave vector Q. The reason is 
that a diverging effective mass (oc 1/Z) causes the Lan- 
dau parameters (dimensionless Fermi-liquid interaction 
components) of all channels except the critical one to di- 
verge as well, thus rendering the corresponding partial 
susceptibilities finite. Acting in the crossed particle-hole 
channel, these divergent interactions cause divergent ver- 
tex corrections oc 1/Z. This means that a fermion-boson 
model in which only critical bosons are included, which 
has been used frequently at weak coupling, is no longer 
sufficient at strong coupling. The interplay of spin fluc- 
tuations and fermionic excitations has also been consid- 
ered in 1/N expansion by Abanov and Chubukov [TU] 
and a renormalization group formulation has been given 
by Metlitski and Sachdev [TT]. However, these authors 
did not consider the quasiparticle renormalization in the 
strong coupling regime. 

In the context of heavy fermion compounds, the quan- 
tum critical point has been often associated with a break- 
down of the Kondo effect, and therefore a breakdown of 
the picture of heavy quasiparticles [T2TTH] . In this sce- 
nario, it is assumed that at the critical point the en- 
ergy scales of the Kondo effect and the exchange interac- 
tion between the localized /-spins (in the absence of the 
Kondo effect) are approximately equal. Some of these 
scenarios have been developed enough to allow compar- 
ison with experimentally observed critical exponents, in 
particular for CeCus.gAuo.i and for YbRh 2 Si2 [15l H6]. 
However, experimentally the quasiparticle mass does not 
appear to be drastically reduced (by orders of magni- 
tude) when the QCP is approached, as would be expected 
if the Kondo effect were to be suppressed. We argue 
here that in the cases we consider the Kondo effect, or 
more precisely the heavy quasiparticle picture, remains 
intact. However, the quasiparticles experience an AFM 
spin-exchange interaction responsible for the ordering of 
their spins in the AFM state. In other words, we propose 
that the ordered state is an itinerant heavy-quasiparticle 
SDW state, at least in the neighborhood of the critical 
point. The resulting small ordered magnetic moment is 
in agreement with observation. Roughly speaking, the 
ordered moments provide a magnetic field acting on the 
Kondo ions. As long as the /-electron Zeeman splitting 
caused by this field is small compared to the Kondo tem- 
perature, the Kondo effect is only weakly suppressed. 

In this report, we present a semi-phenomenological 
theory of the scaling behavior near an AFM QCP within 
an Anderson lattice model. As discussed above, we 



show how energy fluctuations may lead to a momentum- 
independent critical quasiparticle self energy. The feed- 
back of the critical quasiparticle properties (the ^-factor) 
into the spin and energy fluctuation spectrum leads to a 
self-consistent equation for the quasiparticle self energy 
and effective mass m*/m cx Z^ 1 . This allows a strong- 
coupling solution in the form of a fractional power law 
Z{ui) cx ui a with a = 1/4 for 3D and a = 1/8 for 2-D 
AFM spin fluctuations. The free energy obeys scaling 
characterized by fractional power laws. The dynamical 
structure factor is found to satisfy uj /T-scaling within the 
critical cone. Comparison of our theory with experimen- 
tal data shows detailed agreement. 

The above scenario depends sensitively on the detailed 
nature of spin fluctuations in a given system. For ex- 
ample, 3D AFM fluctuations do not lead to true critical 
behavior, i.e. a Gaussian fluctuation theory is applicable 
[U [2] , provided the effective mass enhancement by criti- 
cal fluctuations is not too large. We argue below that in 
CCA there is a wide region of quasi-2D antiferromagnetic 
fluctuations which gives rise to a substantial enhance- 
ment of the effective mass and may drive the system into 
a strong-coupling regime of 2D or 3D antiferromagnetic 
fluctuations. 

Critical quasiparticle picture. Our starting point is 
a heavy Fermi liquid as found for an Anderson lat- 
tice model of correlated /-electrons (on-site interaction 
U) hybridizing with conduction electrons. The energy 
scale of the heavy-fermion band is given by the "coher- 
ence temperature" T co h, which is well above the tem- 
perature regime for which "non-Fermi liquid" behavior 
is observed near a QCP. On top of the heavy-fermion 
quasiparticle renormalization, the critical fluctuations 
cause a further renormalization on which we focus below. 
The single-particle Green's functions may be decomposed 
into a quasiparticle term and an incoherent contribution, 
G(k, uj) — ZG qp + G mc , where the quasiparticle weight 
factor Z is defined by Z" 1 = 1 - dReE(u>)/du. The 
quasiparticle Greens function is given by G qp (k, uj) — 
[uj — Ek — ir] _1 , with Ek = (m/m*)vp(k — hp), where 
vp is the Fermi velocity of the heavy-fermion band and 
the quasiparticle width is T = ZlmX^Ek). 

The condition for the quasiparticle picture to be valid 
is T < \Ek[- In the Fermi-liquid regime, T = c(Ek) 2 <C 
|i£fc| in the limit Ek —> 0. Here, we argue that the 
quasiparticle stability condition may be even satisfied 
in non-Fermi liquid situations. We extend the usual 
quasiparticle picture by recognizing that the parameter 
Z = m/m* depends on the energy scale, Z = Z(ui) — 
1/[1 — dY,(u))/du)]. It is important to observe that the 
(retarded) self energy is an analytic function in the up- 
per half plane, so that the real and imaginary parts of 
any nonanalytic term (in the lower half plane) are locally 
connected. Even in a true non-Fermi-liquid phase with 
E(u>) cx -i^u;!) 1 "", a < 1 so that ImE(w) oc Re£(u;) cx 
H 1 -" and Z cx (\E k [) a , one finds T/\E k \ = tan(f a) < 1 
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for < a < 1/2. In this case, although Z = at the 
Fermi surface, the spectral function for non-zero excita- 
tion energy may be peaked sharply enough to separate a 
quasiparticlc contribution from the incoherent part. 

Spectrum of critical fluctuations. We assume that in 
the paramagnetic phase of a metal close to an antiferro- 
magnetic quantum critical point, the self energy for the 
single-particle Green's function is determined by the in- 
teraction with magnetic fluctuations. Here, we take the 
imaginary part of the renormalized retarded dynamical 
spin susceptibility for wave vectors near the AFM order- 
ing wave vector Q to be of the form 



Imx(q, v) 



(N /Z)(v/v* F Q) 
[r + Zcfe } 2 + {v/v* F QY 



(1) 



where q is measured from the ordering wave vector Q. 
Here, Nq is the bare and Nq/Z the renormalized density 
of states at the Fermi surface, v* F — Zvf is the renormal- 
ized quasiparticle Fermi velocity, £o — hp is the micro- 
scopic AFM correlation length and the control parameter 
r = r + T/Z, where r (H) = 1 + F(Q, H) oc H/H c - 1 
and T/Z is the temperature-dependent shift of the pa- 
rameter r [HGI]. Here, F(Q,H) is a dimensionless gen- 
eralized Landau parameter ,which tends to —1 at the 
critical point. H is representative of the control param- 
eter that tunes to the phase transition, and thus may be 
magnetic field, pressure, doping concentration, etc. The 
factor Z multiplying q 2 in Eq. (1) is a crucial feature of 
our theory, see [3] ■ 

We now define an energy-fluctuation propagator 
Xe{k, v) by combining two spin fluctuation propagators, 
in the form 



ImxE(K,v) = ^2 G fc+9l G fc+(Zl _ K Imx(qi,^i) 



(2) 



xlmx(qi -K,U! -v)[b(vi -v) -b{v x )], 

where b(v) is the Bose function. The Green's functions 
Gk+qn Gk+ qi ~n are off-shell and may be replaced by 
Performing the momentum integration by Fourier trans- 
form, we get 



Z 2 



d/2+l 



1 



(3) 



where we defined ujq ~ vfQ ~ ep and used N 1 rs ijJqS,o- 
Quasiparticle self-energy. Now we set up a self- 
consistent determination of the quasiparticle self energy 
via the leading term in a skeleton graph expansion in 
terms of the boson propagator \e- The imaginary part 
of the self energy is given by 



/j 
— ImG(k + q, u + v) 
n q 

x Imx B (q, v)[b{v) + f(u - v)] 

« W o^VM)Z 2 ) d - 1/2 cx | w |d-V2-«(2d+i) ! 

(4) 

where we used J dcos,9 q: klmGk +q rj (1/Zvpq), the di- 
mensionful (Landau) interaction vertex is given by oc 
(ZNq)^ 1 , and f(oj),b(oj) are Fermi and Bose functions, 
which at low T confine the ^-integration to the inter- 
val [0,w]. In the last equation we assumed the power 
law Z{uj) cc \uj\ a . The scale dependent contribu- 
tion to ReE(w) follows from analyticity as ReE(w) cx 

_| w |e2-l/2-a(2d+l) a %-2 (^^2^-1/2 _ This leadg to 

the self-consistency relation for Z(u>) 



Z^iuS) = l-dReZ{uj)duj w l+Z~ 



2d-l 



d-3/2 



(5) 



We now explore the consequences of the scale depen- 
dent Z. In general, the u> and T dependence of Z is 
obtained by substituting ^/u 2 + cT 2 for ui, where c is 
a constant of order unity. For frequencies less than 
the temperature, we may replace uj by T . As long as 
Z- 2d - 1 (T)(T/u) ) d - 3 / 2 < 1 for any T, the system will 
be in the Gaussian fluctuation regime all the way down 
to the critical point. If, however, the initial value of 
Z _1 (T), when one enters the AFM fluctuation regime, is 
sufficiently large, such that Z- 2d - 1 ( T )( r / w o) d_3/2 > 1, 
a new regime is accessed, which is of a strong-coupling 
nature. We find the characteristics of this new regime 
within the present approximation by solving the self- 
consistent Eq. (5), to get 



Z(T) oc (T/looY 1 - 3 ^^ 2 . 



(6) 



The exponent a is therefore found to be a = 1/4 for 3D 
and a = 1/8 for 2D fluctuations. 

In the case of only AFM fluctuations in a clean system, 
it is difficult to satisfy the strong-coupling condition of 
sufficiently large Z~ l (T) . However, if on the initial ap- 
proach to the critical point, fluctuations dominate that 
lead to a growing Z _1 (T) with decreasing T, the condi- 
tion may be met at some point. The precise crossover 
point is determined by the crossover of these precursor 
fluctuations to the critical AFM fluctuations and by the 
condition above that leads to Eq. (4). As discussed in 
[H [5] impurity scattering helps to enhance the effect of 
AFM fluctuations on Z(T) . In addition, there are clear 
indications in the data on YbRh2Si2 of 3D FM fluctua- 
tions, p] In that case, one finds Z' 1 ^) oc hx(T /T), so 
that Z~ x grows as T — > . 

Role of non- critical fluctuations and vertex corrections. 
An effective mass that diverges uniformly over the Fermi 
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surface i.e. Z — > seems to cause diverging susceptibili- 
ties in all channels, since in Fermi-liquid theory the par- 
tial susceptibilities are given by xi = Nq/{1 + Fi), with 
A^q = Nq/Z the renormalized DOS and F[ the Landau 
interaction parameter in channel I . Anticipating that all 
susceptibilities except the critical one stay finite as the 
critical point is approached, we are forced to conclude 
that the Landau parameters diverge Fi cx 1/Z , com- 
pensating for the divergent DOS. This has been shown 
within microscopic theory in at least one case, a d-wave 
Pomeranchuk instability [19j . The diverging Landau in- 
teraction gives rise to a singular correction of the spin 
vertex at q = Q , Fg cx l/Z , when substituted in the 
crossed p-h channel. The two vertex corrections at both 
ends of the particle-hole bubble Xoil, w ) = GG cancel 
with the qp weight factors of the two Green's functions, 
in analogy to the known behavior at q = . The exis- 
tence of these vertex corrections is an essential ingredient 
of our theory. 

Dynamical scaling. The critical behavior of the spin- 
excitation spectrum as discussed above, see Eq. (1), is 
given by (employing units (u)o,£,q 1 ) for ( u} - ( l) 



Imx(q, w ) (X 



w/Z 2 



[r + Zq*] 2 + (Lj/Zy 



(7) 



By equating the terms Zq 2 and uijZ in the denom- 
inator, with Z oc ui a , the dynamical critical exponent 
is found as z — 4d/3 . The correlation length £(H,T) 
in the Fermi liquid (quantum disordered) regime follows 
from r ~ Zq 2 ~ q 2+az as £(H,T = 0) ~ (H - H c )~ v ', 
where v = 3/(3 + 2c?) . In the critical regime, the lead- 
ing temperature dependence of the control parameter 
r(T) = r a + 5r(T) is generated by the fluctuation interac- 
tion (c/) 4 ) - term in lowest approximation: Sr(T) cx T 1 ~ a . 

It follows that £(ff c ,T) ~ T- 1 ' 2 . The boundary of 
the critical region - the critical cone - in the (H , T) 
phase diagram is found from £(H,T = 0) = £(H C ,T) as 
T ~ \H - H c \ zl/ . Then £(r ,T) has the form f(r ,T) = 
T" 1/z 3(r T- 1 / l/z ). The (u,T) dependence of Z may be 
accounted for with the form Z(w, T) cx [lu 2 + (aT) 2 ] Q / 2 = 
T( z ~ 2 )/ 2z £(w /T). Here, a is a constant (perhaps n). 

Then 



Imx(q,w) cx T 



u/TQ 2 



which shows the following general scaling relation 



(8) 



Im X (q, W ) cxT- 1 ^-,^;^ 1 /^). 



(9) 



Inside the critical cone, where we may set r = and 
cx T 1 / 2 , with x — w£ 2 — lj/T, and y — q^ , we find 
the scaling form 



Imx(q,w) ~ £ z ®(x,y) , 

x(- 2 {x) 



<f>(x,y) 



{l + c q Cy 2 ) 2 + {x/Cr 



(10) 
(11) 



where ((x) = cz[x 2 + a 2 ] Q / 2 and c q , cz are constants. 

Inside the critical cone and at the ordering wave vector 
the spin excitation spectrum obeys w/T-scaling 



Imx(Q,w) 



W /TC 2 



T 1 + {u/TO 2 



(12) 



A comparison of this scaling form with neutron scat- 
tering data on CeCuAu is shown in Fig. 1. 




E/T 



FIG. 1. Inelastic neutron scattering: Comparison of theory 
Eq. (12) and experimental data [20] for CeCu6-ic Au^ at the 
critical concentration x = 0.1 

Free energy. The critical part of the free energy den- 
sity may be derived from the expression for the entropy 
density in terms of the self-energy [2Tj : 



S 
V 



2tt 



JV(0) 



dwuj 



T 2 cosli 



[w-ReE(w)] (13) 



IT 



Substituting the critical self energy found above and in- 
tegrating over temperature we find the scaling form of 
the free energy density in the case of d-dimensional spin 
fluctuations (d = 2, 3) in a 3D metal 



f(H,T)=r {2d+1) ^f(r e / \TC) 



(14) 



The correlation volume V c ~ £( 2rf+1 )is understood as 
follows: The underlying critical degrees of freedom in 
the very low temperature Anderson lattice picture are 
the heavy fermonic quasiparticles described by the prop- 
agator G(k,ui)^ 1 — ui — vk — with cx cj 1_q . 
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Therefore, their dynamical exponent Zf — 1/(1 — a) and 
their dimensionality is dt = 1 [22]. The entropy of the 
system is determined by hyperscaling for the fermions: 
S oc T d f/ Z f. Therefore, since £ oc T _1 / z , where z = 4d/3, 
the free energy density / oc T 1+d f' z J ~ £2d+i^ jj er6j z 
and ci are the bosonic ones discussed above. 
The specific heat coefficient follows as 



C/T oc 



ji(2rf+i 2z)/z^ cr itical regime 



1 , Fermi liquid regime. 

For the critical part of the magnetization we find 



(15) 



M-M(H C ,Q) 



-T, 



«(2d+l)-l 



critical regime 
, Fermi liquid regime. 

(16) 



The susceptibility has a critical part 

( —T 1 !- 1 /^ critical regime 

X- X{H C ,0) Ct < „(2d+l)-2 „ . r 

I — r , lermi liquid regime. 

(17) 

Using <9Af /<9T = —const, we find the Grueneisen ratio 
in the critical regime: 

r c = -^«£«T<-*-iv (18) 

while in the Fermi liquid regime, we have the universal 
result 



Tg = -h%- c > G r = -(z- 2d ± 3 )u (19) 

The critical scaling of the transport properties is ob- 
tained by observing that the quasiparticle relaxation rate 
scales as T oc £~ z $r{rot; 1 ^ , T£ z ). Then the resistivity 
behaves as 



p-p(H c ,0) oc — r 
m 



rp(z+2)/2z 
_| ro |(2-3^)t//2y2 



(20) 

critical regime 
Fermi liquid regime. 



For the thermopower, S we get, using S oc m*T, S oc C 
in both the critical regime and the Fermi liquid regime. 

Comparison with experiment. The above results, in 
the case of 3D fluctuations, are identical to the ones ob- 
tained by two of us previously [5] [H]> where impurity 
scattering was invoked to give rise to a critical, weakly 
momentum dependent self energy. In particular, the Z- 
exponent a = 1/4 found there, and the ensuing critical 
indices z = 4, v = 1/3 are the same as the ones found 
in the present work for the clean system. The excellent 
agreement of the theory 8 , 9 with the experimental data 



on YbRli2Si2 in the regime close to the QCP therefore 
applies to the present theory as well. In contrast to that 
earlier work, the present results do not depend on the im- 
purity concentration. Indeed, in experiment, the critical 
parts of e.g. the specific heat do not show a dependence 
on the impurity content of the sample. 

We turn to a different case, CeCug-zAoj, for which 
a QCP has been found at the concentration x = 0.1 at 
ambient pressure in the absence of a magnetic field and 
at slightly different concentrations at a critical pressure 
and/or an applied critical field. As suggested by the neu- 
tron scattering data, the magnetic fluctuations there ap- 
pear to be quasi-two-dimensional. We therefore compare 
our results for d = 2 with the available data. In Fig. 
1, we have already compared theory and experiment for 
the dynamical spin susceptibility. In Fig. 2, we show the 
specific heat data [53] in comparison with C/T oc T -1 / 8 
as obtained above in Eq. (15) 



• data 
C/T = -3 + 4T" 1 / 8 




T{K) 

FIG. 2. Specific heat: Comparison of theory Eq. (15) and 
experimental data [23] for CeCue-^Aua; at the critical con- 
centration x — 0.1 

The resistivity result p{T) ~ p(0) oc T 7 / 8 is fitted to 
the data in Fig. 3 

Our prediction for the uniform magnetization is 
M(T) = M(Q) - aT + bT 2 , from Eq. (16) augmented 
by a Fermi liquid correction oc T 2 . It fits the data [23] 
well, as shown in Fig. 4. 



Conclusion. 

We presented a semi-phenomenological theory of quan- 
tum criticality near a critical point that separates antifer- 
romagnetic and paramagnetic phases of a metal. Start- 
ing from the assumption that the Landau quasiparticle 
effective mass diverges on approaching the critical point, 
giving rise to critical quasiparticles, we identified the cor- 
responding renormalization of the dynamical spin suscep- 
tibility near the antiferromagnetic wave vector. Critical 
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FIG. 3. Resistivity: Comparison of theory Eq. (20) and ex- 
perimental data [24] for CeCue-xAu^ at the critical concen- 
tration x = 0.1 
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FIG. 4. Uniform magnetization: Comparison of theory Eq. 
(16) and experimental data [23] for CeCue-xAu^ at the crit- 
ical concentration x = 0.1 



contributions to the electron self energy induced by an- 
tiferromagnetic fluctuations in a clean system are known 
to be strongly anisotropic (confined to the "hot spots"). 
However, impurity scattering may be shown [5] [5] to dis- 
tribute the effects of critical scattering all over the Fermi 
surface. Here we have shown that even in a clean sys- 
tem, the critical antiferromagnetic fluctuations give rise 
to a critical self energy uniformly over the Fermi surface. 
This is because the magnetic fluctuations generate en- 
ergy fluctuations, which diverge in the long- wavelength 



limit. The scattering of quasiparticles off these energy 
fluctuations leads to a contribution to the effective mass 
which is nominally proportional to a positive power of 
energy, but is strongly enhanced by factors of the effec- 
tive mass itself. This leads to a self-consistency relation 
for the effective mass, the quasiparticle Z-factor, which 
may have a strong-coupling solution provided the initial 
value of Z at the appropriate high-energy scale is suffi- 
ciently small (as may be caused by fluctuations leading 
to a weakly-diverging effective mass). While the criti- 
cal quasiparticles living on the "cold" parts of the Fermi 
surface dominate most of the observable quantities, the 
"hot" quasiparticles may be shown to be even more singu- 
lar, e.g. in d = 3, we find the equivalent of the exponent 
a to be =1/2. In this context, it is interesting to observe 
that the observed critical behavior of CCA depends on 
whether the QCP is tuned by varying the Au concentra- 
tion, the pressure, or the magnetic field. It is conceiv- 
able that in the case of magnetic field tuning the precur- 
sor fluctuations necessary to access the strong coupling 
regime are too weak so that the system remains in the 
weak-coupling regime, as is apparently observed. A fur- 
ther condition for the applicability of the self-consistent 
solution is that the effective dimension of the bosonic 
fluctuations, z + d, is sufficiently above the upper critical 
dimension of the appropriate field theory (e.g. </> 4 the- 
ory), such that boson-boson-interaction effects may be 
neglected. 

Application of our theory to the cases of 3D and 2D 
fluctuations in a three-dimensional metal leads to a crit- 
ically diverging effective mass m* cx T~ a with a = 1/4 
(3D) and a = 1/8 (2D), in good agreement with experi- 
mental data on the two heavy-fermion compounds YRS 
and CCA, where neutron scattering showed the presence 
of 3D and 2D AFM fluctuations. Our theory obeys hy- 
perscaling, taking into account the scaling of the critical 
fermions. Further comparisons of our theory with exper- 
imental data on YRS and CCA show good agreement. In 
particular the universal behavior of the Griineisen ratio 
in the quantum-disordered regime measured in YRS is in 
excellent agreement with our result [5]. 

We also note that we have found that the correlation 
length inside the critical cone grows unusually slowly 
with decreasing temperature, which appears to have been 
observed in neutron scattering and has sometimes been 
interpreted as supporting the notion of a "local quantum 
critical point." 

The results of this paper can, alternatively, be derived 
in terms of a renormalization group (RG) approach [25j . 
The RG analysis demonstrates that the theory is for- 
mally an expansion around the lower critical dimension 
where a marginal Fermi liquid behavior occurs and the 
notion of singular quasiparticles is well-defined. In our 
case, this corresponds to did = 3/2, where the exponent 
a vanishes. The physics discussed in this paper then ap- 
pears as an unstable fixed point. It implies that in the 
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limit where Z — > 0, the perturbation theory used here is 
ultimately expected to break down. Nevertheless, we ex- 
pect the power laws obtained here to govern the relevant 
physics for a wide range of energies and temperatures. 

Finally, we emphasize that we assumed the heavy 
quasiparticles to be robust, though modified by scatter- 
ing from critical spin fluctuations. There is no break- 
down of the Kondo effect nor an associated collapse of 
part of the Fermi surface in our scenario. Experimental 
features, such as the crossover behavior observed in trans- 
port properties across the "T*-line" in the T — H phase 
diagram of YRS, may be accounted for by a change in 
quasiparticle scattering strength associated with thermal 
activation of the (ESR) spin resonance [26] . 

While the good agreement of our theory with exper- 
iment across the board is encouraging, there is a need 
to put the phenomcnological assumptions it involves on 
a firm microscopic basis. Work in this direction is in 
progress. 
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